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-Computer program for centrifugal distortion in asymmetric 


3 _ top rotational spectra 


4 _ By Gunnar ERLANDSSON 


The electronic computer Besk has previously been applied to the analysis of 
asymmetric top rotational spectra in the rigid case (1) (2). In the present work 
the use of Besk has been extended to a first order treatment of centrifugal dis- 
tortion in such spectra. The program has been tried on a group of Q-branch 


lines in the spectrum of formaldehyde and has given a good representation of 


observed distortions. The results obtained are consistent with results obtained 
previously by other methods. 


Method of calculations 


As shown e.g. in (3) the influence of centrifugal distortion on the rotational 
energies may be represented to first order of approximation by an expression of 
the form 
W = Wot D, Wo+ Dod (J +1) Wot D3 J? (J +1)? + Did (J +1) < P25 + 

+.Ds<Pa> + De Wo<Pa> [1] 


where W, is the rigid rotor frequency, P, is the component of angular momentum 


along one principal axis of the molecule, and the brackets <> denote average 


_ values for the energy level in question. 


The distorted value of a line frequency is obtained as the difference between 
two equations of the form [1]. If a sufficient number of rotational lines are 
known, this formula gives a system of equations from which the constants D, 
to D, may be determined. For this purpose, however, one must know the values 
of W,, <P> and <P%> for all levels involved. These values may all be calculated 
from known values of the rotational constants A, B and C' of the molecule. 

The general procedure for the application of formula [1] is therefore to start 
from values of A, B and C, obtained from reasonably undistorted lines or from 
suitable approximation methods, and to calculate Wo, <P%>, and <P%> for all 
levels of interest. The differences A W,=%, A Wo etc. are then formed, giving a 
system of equations for the D;, which is finally solved by the least squares method. 

The present program consists of three parts. Part 1 calculates the reduced 
energies, H(b), together with <P?) and <P%> for ali levels of each sub-matrix 
involved in the transitions considered. Part 2 calculates the complete energies 
W, for the particular levels involved, and sets up the resulting system of equations. 
Part 3, finally, gives the solution of this system. 
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Details of the program 


Part 1 of the program sets up the elements of the energy sub-matrix con- 
sidered, using the asymmetry parameter b= (C—B)/(2A BG) and calculates 
the reduced energies #(b) by means of the continued fraction method (4). The 
matrix representation employed is diagonal in the prolate symmetric rotor case 
(6=0). The calculations of #(b) start from approximate input values and proceed. 
to:'a relative accuracy better than one part in 104. 

A difficulty often encountered in calculating reduced energies by the continued. 
fraction method is sporadic divergence of the fractions for special values of the 
asymmetry parameter. This situation may always be remedied by choosing a. 
development starting with another diagonal element in the matrix. The present. 
program considers all developments for each case simultaneously, and selects 
the one giving the best convergence. This procedure results in a comparatively 
simple and effective program avoiding the administration procedures for changing 
the development form. Changes in the approximations tried are at first limited. 
to absolute values smaller than unity. Refinements are then calculated by the 
“regula-falsi’? method. 

After evaluating all H(b):s for one sub-matrix and storing the results in the 
memory, the program proceeds to calculate the elements of the transformation 
matrix 7’, which diagonalises the energy matrix. The transformation 7 is then 
applied to the diagonal matrices K* and K*, which represent P% and P* in the 
original system. The diagonal elements of the matrices thus transformed give the 
desired averages <PZ> and <P{>. These are printed in an output table together 
with corresponding values of H(b). If desired the program then proceeds auto- 
matically with other matrices or other b-values. 

The 7-matrix is evaluated according to the procedure given by King, Hainer, 
and Cross (4). This requires a very high degree of accuracy in the values of the 
reduced energies and in the arithmetic operations involved. A small error in the 
first elements of a column may grow to considerable magnitude at the bottom 
of the column by successive multiplications with elements in the order of ag 
where K takes on alternate integer values up to J. These difficulties are espe- 
cially pronounced near the symmetric case b=0. 

In the present program the accuracy obtainable is limited by the word-length 
of the machine, which is 40 binary digits. Calculations are performed in a float- 
ing-point manner conserving the relative accuracy as far as possible. For future 
work an improved program using the double word-length is being planned. 

In the case of formaldehyde, which is a nearly symmetric rotor, the straight- 
forward calculation of 7 could not be used for J>13. For higher J:s errors 
occurred in <P%> and <P%), which could be observed by checking the traces of 
the matrices. The errors in the 7'-elements occurred at first far below the main 
diagonal. An improvement was obtained by setting such elements equal to zero, 
as soon as they could be estimated to be unimportant. This procedure gave 
satisfactory results up to J=17 for the formalhyde case. 

Part 2 is a straight-forward arithmetic and administrative program evaluating 
the coefficients of the, desired system of equations. Part 3 is a standard program 


for the solution in Besk of over-determined systems of linear equations by the 
least squares method. 
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Table 1 


a ES oe, irked by 
Observed | Rigid rotor Centrifugal | Calculated— 


Transition frequency frequency distortion distortion 
* Me Me Mec Me 
ee ee 1 a eee 
21.2-21.1 14488.65 14488,22 — 0.43 — 0.56 
31.3-31.2 28974.85 28972.17 — 2.68 — 2.62 
41.4-41.3 48284.60 48276.76 — 7.84 ma O9 
51.5-51.4 72409.35 72391.71 — 17.64 Sta 
72.6-72.5 8884.87 8871.95 — 12.92 = 13.19 
82.7-82.6 14726.74 14702.20 — 24.54 — 24.70 
92.3—92.7 22965.71 22922.12 — 43.59 — 43.39 
14342-14341 7892.03 7834.86 — 57.17 — 57.04 
153.13-153.12 11753.13 11662.68 — 90.45 — 90.43 
163.14-163.13 17027.60 16888.33 — 139.27 — 139.34 
173.15-173.14 24068.31 23859.09 — 209.22 — 209.22 


Practical results 


The program has been tried on a group of eleven branch lines in the micro- 
wave spectrum of formaldehyde. These lines were reported by Lawrance and 


_ Strandberg (5), who accounted for the centrifugal distortion by means of a semi- 


classical method. The distortions involved are comparatively large, ranging from 
0.4 Mc for J=2 to 209 Mc for J=17. Lawrance and Strandberg also considered 
six lines of higher J-values, but for reasons explained above these could not be 
included in the present treatment. 

The results obtained for the individual lines are shown in Table 1, which gives 
observed frequencies (y), calculated rigid rotor frequencies (y)), centrifugal distor- 
tions (v—v), and distortions calculated from eq. 1 with the constants obtained 
as results from Besk. 

The rotational constants used for the rigid rotor case are taken from (5) and 
(6) and are: 


A = 282.106 kMc 
B= 38.83407 kMc 
C= 34,00437 kMc 


The distortion constants for eq. 1 derived by Besk are: 


D, = — 0.28627 - 10-8 (+£0.039 - 10-8) 
D,=  0.12209- 10-4 (£0.028 - 10-4) 
D,= —0.46625 - 10-3 (+£0.081 - 10-3) 
D,=  0.14253-10-4 (+0.047- 10-2) 
D,= 0.15399 - 10-3 (+0.022 - 10-8) 


The constant D, (corresponding to D,;) drops out from the equations for the 
Q-branch lines. The uncertainties given in parenthesis are those pertaining to the 
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"From the Dys 
Ppecdived to the 54 +1)! 
tween the various sets of co 


This may be compared to the value 
symmetric top approximation to the R. 
very high in neither case, the agreement m 

As mentioned above an improved Rie € 
later on, and it is hoped that this will make} ible he 
the lines observed. It has therefore not been considered wort 
further with the analysis at present. Otherwise, some minor a 
rigid rotor rotational constants should be tried, and the procedure be 
for a best fit. s 
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